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We investigate the effects of the underlying quantum geometry in loop quantum cosmology 
on spacetime curvature invariants and the extendibility of geodesies in the Bianchi-I model 
for matter with a vanishing anisotropic stress. Using the effective Hamiltonian approach, we 
find that even though quantum geometric effects bound the energy density and expansion 
and shear scalars, divergences of curvature invariants are potentially possible under special 
conditions. However, as in the isotropic models in LQC, these do not necessarily imply a 
^vq , physical singularity. Analysis of geodesies and strength of such singular events, point towards 

a general resolution of all known types of strong singularities. We illustrate these results for 
O4 the case of a perfect fluid with an arbitrary finite equation of state w > — 1, and show that 

■ curvature invariants turn out to be bounded, leading to the absence of strong singularities. 

Unlike classical theory, geodesic evolution does not break down. We also discuss possible 
generalizations of sudden singularities which may arise at a non-vanishing volume, causing 
a divergence in curvature invariants. Such finite volume singularities are shown to be weak 
and harmless. 
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I. INTRODUCTION 

(N; 

Theorems of Penrose, Hawking and Geroch prove that existence of space-like singularities is a 
\ generic feature of Einstein's theory of general relativity (GR) These events can be characterized 

in various ways, such as in terms of divergences of curvature invariants, breakdown of geodesies, 
tidal forces becoming infinite etc. They signal that the limit of applicability of classical gravity 
has been reached. It is believed that new physics from a quantum theory of gravity will provide 
insights on the resolution of singularities, leading to a non-singularity theorem. However, before 
we implement the techniques of quantum gravity to attempt the resolution of singularities, it is 
important to distinguish which of these are harmful, and which are harmless. Neither a divergence 
of curvature invariants implies an end to geodesic extendibility nor geodesic incompleteness nec- 
essarily implies existence of a physical singularity. To understand the nature of a singularity one 
needs to analyze its different properties in unison. In particular, it becomes important to examine 
the strength of singular events - whether they lead to an inevitable destruction of in-falling objects 
or allow passage of sufficiently strong objects. The former characterize strong singularities, and 
the latter, weak singularities [244] ■ It has been conjectured that strong singularities correspond to 
events which can not be geodesically extended 0, @] , and thus represent the true singularities that 
need to be resolved. 

In recent years, methods of loop quantum gravity (LQG) have been successfully used in the 
quantization of homogeneous spacetimes to address the issue of singularity resolution @, 0] • In 
this canonical framework, known as loop quantum cosmology, a complete quantization of various 
isotropic spacetimes has been performed for several matter models [9H19[]. These investigations 
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show that in various loop quantized spacetimes, the big bang singularity of GR is replaced by 
a big bounce which occurs when spacetime curvature approaches Planck regime. 1 These results 
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1 The existence of bounce has recently been demonstrated in the consistent probabilities framework applied to LQC 



have been extended to the anisotropic spacetimes where the quantum evolution has been shown to 
be non-singular [22-25], and rich physical implications have been derived 2^-31]. The resulting 
physics of LQC stands in striking contrast to those in Wheeler-DeWitt quantum cosmological 
models where singularity resolution has remained a formidable challenge to overcome. A key aspect 
of LQC which separates it from previous attempts is the underlying quantum geometry inherited 
from LQG. It has various profound implications. Not only it is responsible for the quantum 
Hamiltonian constraint to be a quantum difference equation in the geometric representation in 



LQC, it also leads to upper bounds on energy density [ll|, |l5|, [22j, |3l|, |32| and expansion and shear 



scalars of the congruence of geodesies [271 . 
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Existence of these bounds provide important 
insights on singularity resolution and also the new physics in loop quantum spacetimes. As an 
example, a spatially flat isotropic universe sourced with a matter content with equation of state 
satisfying weak energy condition (WEC) bounces in LQC when the upper bound for energy density 



is saturated 11 



34i | . Thus, avoiding the big bang/crunch singularity. On the other hand, if it is 
filled with matter which violates WEC,it recollapses when the maximum of energy density is 
reached, and avoids big rip singularity 0,113. 



An important question is whether the singularity resolution in LQC is a generic phenomena. In 
particular, whether all strong singularities of the classical theory are resolved in LQC. This ques- 
tion was recently addressed by the author in the spatially flat isotropic model using the effective 
dynamics of LQC resulting from an effective Hamiltonian [33]. The effective Hamiltonian approach, 
based on the geometrical formulation of quantum mechanics [s^], provides an invaluable tool to 
understand the underlying discrete quantum evolution in terms of variables in the continuum ef- 
fective spacetime. The resulting physics obtained from the effective dynamical equations has been 
extensively tested to confirm with the underlying quantum difference equation for isotropic and 
anisotropic models 11, 12, 14, 19, 39, [13]. The modified dynamics from the effective Hamiltonian 



has been used to prove that all strong singularities of the spatially flat isotropic LQC are resolved. 
These include big bangs/crunches and also finite volume singularities, such as the big freeze singu- 
larity [4l|, which can arise in models with a generalized equation of state. The effective spacetime 
of spatially flat isotropic LQC turns out to be geodesically complete [33)] • Investigations carried out 
for spatially curved isotropic models in LQC, strongly suggests similar results [42J. It is interesting 
to note that these studies also find events where curvature invariants diverge in isotropic LQC. 
Such events correspond to the sudden singularities, which occur at a finite volume and energy 
density but have a divergent pressure and have been recently studied in GR 4314461] 3 . All of the 
events where curvature invariants diverge in isotropic LQC, turn out to be weak singularities which 
are harmless. 4 

The goal of this work is to analyze these questions for Bianchi-I models in LQC using the 
effective Hamiltonian approach. Due to an interplay of the non-vanishing Weyl curvature with the 
Ricci curvature, Bianchi-I spacetimes provide one of the simplest expositions of the rich structure 
of singularities in GR, and thus an opportunity for a non-trivial and important generalization of the 
results on generic resolution of strong singularities so far proved in isotropic models in LQC. For 
perfect fluids, depending on the equation of state of matter and the initial anisotropy, singularities 
can be of the form of a barrel (one of the scale factors is finite and other two vanish), a cigar (one 



[20(| . Similar methods for Wheeler-DeWitt quantization show that the probability for the singularity to occur is 
unity JiJ]. 

The big rip singularity in GR occurs when volume diverges in finite time with a divergences in energy density. For 
an example of a scalar field model leading to such a singularity, see Ref. ^35 
For a discussion of such singularities in modified gravity scenarios, see Ref. 
The first example of such a 'singular' event in LQC was first reported for a scalar field model in Ref. |48|]. It 
was then shown that such events can arise quite commonly in isotropic LQC with a suitable choice of generalized 
equation of state [3^] . 



4711 and reference there in. 
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of the scale factors diverges and other two vanish), a pancake (one of the scale factors vanishes 
and other two approach finite non-zero values) or a point which are isotropic singularities (all 
scale factors vanish) 4^-54]. In contrast to the isotropic models, where for a fixed equation of 
state, behavior of energy density is sufficient to capture the details of the spacetime curvature, in 
Bianchi-I model one must take in to account the important role played by the shear scalar. Roughly 
speaking, it is both the energy density and the shear scalar which determine whether the expansion 
scalar behaves more in an isotropic or in an anisotropic way. This is also reflected in the effective 
dynamics of Bianchi-I model in LQC, where bounces do not occur at maximum value of energy 
density, but when at least the energy density or the shear scalar approach values in the Planck 
regime 27|, 1 3 11 ] - Due to these rich features, analysis of various aspects of singularity resolution, 



such as the strength of singularities in Bianchi-I model using Tipler and Krolak's conditions is 
more subtle than in the isotropic models Qj. To simplify the analysis, we restrict ourselves to the 
study of matter with a vanishing anisotropic stress, and aim to answer following questions related 
to resolution of singularities in the effective dynamics of Bianchi-I LQC. First, are the curvature 
invariants bounded in LQC or as in the isotropic spacetimes they can potentially diverge? If so, 
under what conditions such divergences occur? Secondly, under what circumstances the geodesic 
equations in the effective spacetime breakdown? In particular, can geodesies be extended beyond 
the events where curvature invariants diverge? Finally, do these events correspond to strong or 
weak singularities? 

The analysis carried out in this work and the main results are organized as follows. In Sec. 
II, for completeness, we summarize the dynamical equations obtained from the connection-triad 
variables for the orthogonal Bianchi-I model in GR as well as in effective spacetime description 
of LQC. Here we demonstrate the way universal bounds on energy density, expansion and shear 
scalars of geodesic congruences are reached. (For further details on this part, we refer the reader 



to Refs. [2g . 1271 . l3ll . |32j). These dynamical equations are used to obtain the expressions for 
curvature invariants - Ricci scalar, Kretschmann scalar and the square of the Weyl curvature in 
Sec III. We then find the conditions under which curvature invariants can potentially diverge in 
the effective dynamics of LQC. In Sec. IV, we first analyze the behavior of null geodesies and 
obtain the conditions for which they can break down in the effective spacetime of Bianchi-I model 
in LQC. These conditions show that events where curvature invariants may diverge in LQC are 
not necessarily geodesically inextendible. We then perform an analysis of the strength of singular 
events. Interestingly, the conditions where singularities can potentially be strong turn out to 
be identical to those which determine the inextendibility of geodesies. In Sec. IVC, we discuss 
the physical implications of these results for a (non-viscous) perfect fluid with an arbitrary finite 
equation of state w > — 1. For this model, we show that unlike the classical theory, there exist no 
curvature divergent events in LQC. The effective spacetime in LQC turns out to be geodesically 
complete and devoid of any strong singularities. This is followed by a discussion of the fate of 
certain types of exotic singularities which may arise with a generalized equation of state. In 
particular, we analyze the singularities which may occur at finite volume with a divergence in 
curvature invariants. Such events, which generalize sudden 43J] and big freeze singularities 41[ of 
isotropic models to anisotropic spacetimes, are shown be geodesically extendible weak singularities. 
We summarize the results in Sec. V. 
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II. CLASSICAL AND EFFECTIVE HAMILTONIAN IN ASHTEKAR VARIABLES: 

SOME KEY FEATURES 



We consider a homogeneous (orthogonal) Bianchi-I spacetime with a spatial manifold M 3 . In 
order to introduce a symplectic structure, one needs to introduce a fiducial cell V on the manifold. 5 
The cell has a fiducial volume V = hhhi where l^s denote the coordinate lengths, with respect to 
the fiducial metric q a b on M 3 . The physical spacetime metric for the orthogonal Bianchi-I spacetime 
is given by, 

ds 2 = -N 2 dt 2 + a\ dx 2 + a\dy 2 + a\dz 2 (2.1) 

where TV" is the lapse function. The basic gravitational variables in LQG are the matrix valued 
Ashtekar-Barbero connection A l a and triad Ef. Due to underlying symmetries of the spacetime, 



these are symmetry reduced to connection a and triads pi (where i = 1,2,3) [22j, [56|, 157J]. The 



triads are kinematically related to the directional scale factors as to the triad components pi as 

M = hh a 2 a 3 , \p 2 \ = hh aia 3 , \p 3 \ = hh a 2 a 3 . (2.2) 

The modulus sign arises because of the orientations of the triad. In the following, we choose the 
orientation to be positive and fiducial lengths k to be unity, without any loss of generality. In the 
phase space, the connection and the triad variables satisfy: 

{ci,pj} = SirGjSij (2.3) 

where 7 = 0.2375 is the Barbero-Immirzi parameter in LQG. After imposition of symmetries, the 
only constraint we need to solve is the Hamiltonian constraint. This is expressed in terms of the 
holonomies of connection and triads, and quantized to obtain the physical solutions in LQC. In 
the following, we first discuss the dynamical equations resulting from the Hamiltonian constraint 
in GR for the Bianchi-I spacetime. This is followed by the analysis with the effective Hamiltonian 
of LQC in Bianchi-I spacetime. 



A. Classical dynamics 

For lapse N = 1, the classical Hamiltonian constraint in the symmetry reduced variables q's 
and p^s can be written as, 

N 

U c \ = ~ ^Q^2y ^ ClPl ° 2P2 + C3 ^ 3 ° lPl + C2P2 C3P3 ^ + ^ matt ' ( 2 - 4 ) 

where H m att denotes the matter part of the Hamiltonian constraint. Using, T~L C \, dynamical equa- 
tions can be obtained by solving for the Hamilton's equations for gravitational and matter phase 
space variables (such as 4>,p^, if we consider a scalar field model). In the gravitational sector, these 
equations are obtained as follows: 6 

cfH 

p l = {Pi, n cl } = -8vrG 7 -j^ (2.5) 



5 



The cell plays the role of the infra-red regulator which is removed by taking the limit V — > R 3 . Physical predictions 
must not depend on the choice of this cell. This holds for the analysis in this work, as is for several works in LQC. 
For a discussion on ramifications of these limits and invalidity of unphysical quantizations, see ref. [HH ]. 
The calculation for matter part follows the same strategy, and is simpler. For a scalar field, the matter phase space 
variables satisfy {4>,P4>} = 1. In this case, Hamilton's equations yield Klein-Gordon equation [27J. 
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and 

£i = {ci,H cl } = 8irG 7 -^ , (2.6) 

OPi 

where a 'dot' denotes time derivative with respect to proper time t. Using eq. (12.5p . a relation 
between connection components and the directional Hubble rates, defined as Hi = en /en, follows 

Ci = jkdi = 7-E/jOj . (2.7) 

Considering matter with a vanishing anisotropic stress, i.e. when p(pi,p 2 ,P3) = p(piP2P3), {c%Pi — 



CjPj) turn out to be constant of motion [261. |27|. Using above Hamilton's equations, this in turn 
implies that, 

CiPi - CjPj = V{Hi - Hj) = jKij , (2.8) 

where Kij is a constant anti-symmetric matrix. This feature of GR, that V(Hi — Hj) are constants 
of motion, results in the characteristic dependence of shear scalar on 1/U 2 . To understand this, let 
us recall that the covariant derivative of the unit fluid velocity v a tangent to time-like geodesies, 
can be written as 



J fi;u — v fiu 



c«V (2.9) 



where 9 a p is the expansion tensor and oj a p is the vorticity tensor. 7 Since the considered fluid 
velocity is orthogonal to the spatial hypersurface, vanishes. The expansion tensor can be 
expressed in terms of traceless and trace parts: 

<V = ^(sv + v n vU ) + °> • ( 2 - 10 ) 

where 9 is the expansion scalar: 

9 = V/V = {H 1 + H 2 + H 3 ), (2.11) 
and a^y is the shear tensor, satisfying a^v u = 0. Its magnitude defines the shear scalar a 2 : 

a 2 = ± {(Ht - H 2 f + (H 2 - H 3 ) 2 + (H 3 - H^ 2 ) . (2.12) 

Another measure of anisotropic shear which is used commonly is the scalar S 2 , defined as E 2 := 
cr 2 U 2 /6. Using eq. (|2.8p . it turns out to be 

S 2 = k\ 2 + K23 + • (2-13) 

The scalar S 2 is thus a constant of motion in GR for matter with a vanishing anisotropic stress, 
and thus a 2 oc 1/U 2 in GR. In the classical theory, the expansion scalar 9 (|2.1ip and shear scalar 
a 2 are related to the energy density, defined as p = 'H ma ,tt/V, as 

° 2 8ttG 2 

-9=— p + °- (2 - 14) 

Since, mean Hubble rate H = 9/3, above equation provides a generalization of the isotropic 
Friedmann equation in the Bianchi-I spacetime. Solutions of the dynamical equations have been 



The expansion tensor can be introduced in a similar way for null geodesies [58 
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extensively studied for different matter content 4914521] . These exhibit singular behavior at van- 
ishing physical volume where geodesies break down All the singularities studied so far in the 
classical Bianchi-I model are strong in nature which cause inevitable destruction of all the detec- 
tors falling in them 0, 0]. Late time evolution of the universe, whether it becomes isotropic or 
remains anisotropic depends on the matter content. From (|2.14p . we find that the expansion scalar 
isotropizes, at large volume, if the universe is filled with matter such as dust or radiation or when 
ever the energy density of the matter decays slower than the anisotropic shear scalar a 2 in GR. 



B. Effective dynamics 

Loop quantization of Bianchi-I spacetimes has been rigorously performed for a massless scalar 
field and implications for the quantum theory have been studied recently [2~3 ] (see also Ref. [iil]). 8 
To quantize, the classical Hamiltonian constraint is expressed in terms of the elementary variables 
for the quantum theory: the holonomies of the connection A l a taken over closed loops and the 
fluxes of the triads E% (which turn out to be proportional to triads). The procedure leads to 
a non-local nature of the field strength of the connection in the quantum theory, resulting in a 
quantum difference equation in geometric representation. The quantum difference equation turns 
out to be non-singular and results in a continuous differential equation for Wheeler-De Witt theory 
when spacetime curvature becomes small. An important technique to extract new physics result- 
ing from loop quantization of cosmological spacetimes is the effective Hamiltonian method. The 
effective spacetime description of LQC is derived using techniques of the geometric formulation of 
quantum mechanics ^^]. In this formulation, one treats the space of quantum states as an infinite 
dimensional quantum phase space and seeks a faithful embedding of the classical phase space in 
the latter. The underlying procedure, requires a judicious choice of coherent states, and results in 
an effective Hamiltonian up to well controlled approximations. 9 Using the effective Hamiltonian, 



modified dynamical equations incorporating quantum geometric effects can be derived 6114631] . 
Resulting effective equations have been tested in different models using extensive numerical simu- 
lations, and have been demonstrated to capture various details of underlying quantum evolution to 
an excellent accuracy for isotropic models EMI 

[lil . QjJ and also anisotropic models [jj^, |4(J • In 
the Bianchi-I model, the effective Hamiltonian has been used to understand various physical im- 
plications of the loop quantization 32] and have also been studied in the context of Gowdy 



spacetimes jflij]. As in these works, we would assume that the effective spacetime description for 
the Bianchi-I model to be valid in our analysis. 

For lapse N = I, the effective Hamiltonian constraint for the Bianchi-I model in LQC is given 

by 

1 /sin(/2ici) sin(/Z 2 c 2 ) \ , . 

^eflf = - „ „ 9T , z z P1P2 + cyclic terms + 7i ma tt (2-15) 

8irGj z V \ [i\ fj,2 J 

where jii are proportional to the length of the edges of the loop over which holonomies are evaluated 
and are given by, 



H\ = \J—^-, i«2 = AW — , and /i 3 = Aw— ^- . (2.16) 

V PlP3 V V\P2 



8 For earlier works in LQC on Bianchi-I spacetime, see Refs. [39l. l4p|. 1571. 159. l60l| . 

9 There also exist terms proportional to quantum fluctuations resulting from the quantum properties of the state. 
These turn out be negligible for the phy sical universes which grow to a macroscopic size, as is demonstrated by 
several numerical simulations [1^, [II Gil. Our analysis will assume that such terms can be neglected for the 
effective Hamiltonian for Bianchi-I model in LQC. 
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Here A denotes the square root of the minimum eigenvalue of the area operator in LQG: A = 
4^/3^7^, which is the minimum area to which the loop can be shrunk in the quantum theory 



65[. The matter part in the effective Hamiltonian is treated as obtained from Fock quantization. 10 



The effective Hamiltonian leads to the modified dynamical equations, which can be obtained using 
eqs. (|2.5p and (|2,6p . An immediate consequence of the equation of motions for the triads is that 
the connection components, unlike in the classical theory, are not proportional to dj. Due to this 
reason, the directional Hubble rate for a particular scale factor is determined by all the connection 
components. As an example, the relation for H\ is given by 



Hi 



1 



2 7 A 



(sin (ma - ma) + sin (ma + M3C3) + sin (mci - /X3C3)) 



(and similarly for H2 and #3). Here we have used 



Pi 

Pi = -j (sm(fljCj) + sin( / u fc c fc )) cosQZjCj) 
and the relation between the directional Hubble rate and the time derivatives of triads: 



Hi 



Pj Pk Pi 



J>j Pk Pi , 

with i,j and k taking different values corresponding to the anisotropic directions. 
Though, Cj / jdiHi in LQC, using the Hamilton's equation for connection, 
1 



(2.17) 



(2.18) 



(2.19) 



ci 



2p i7 A 



C2P2 cos(// 2 c 2 )(sin(//iCi) + sin(/z 3 c 3 )) + c 3 p 3 cos(// 3 c 3 )(sin(/xiCi) + sin(^ 2 c 2 )) 



- cipicos(/iici)(sin( / u 2 c 2 ) + sin(/i 3 c 3 )) - mPiPi, 
+ sin(/2ici) sin(/2 2 c 2 ) + sin(^ 3 c 3 ) sin(/2ici 



sin(/i 2 c 2 )sin(/i3C 3 ) 



8vrG7 



P2P3 I P 
Pi 



+ Pi 



dp 
dpi 



and similarly for c 2 and C3, it is straightforward to verify that 

d 



At 



(ciPi - cjpj) = . 



(2.20) 



(2.21) 



Thus, as in GR, (c^Pi — CjPj) are constants of motion in LQC, equal to 7^ (eq. (|2.8p ) |26l . |27I |. 
However, CiPi — CjPj / V(Hi — Hj). Thus, shear scalar a 2 in LQC is not proportional to 1/V 2 . 
In other words, if one defines S^q^ = cr 2 V 2 in LQC, Slq C is not a constant. In the limit when 
spacetime curvature becomes small, S^q C approaches the constant value E 2 of the classical theory. 

Novel properties of the above effective dynamical equations have been studied in various works 
26 29, [3l] . Here we summarize the results for the expansion scalar, energy density and the shear 



In principle, there can be quantum geometric modifications to the matter Hamiltonian, originating from the inverse 
triad terms in the loop quantization. However, these can only be meaningfully defined for compact topologies and 
are significant only at the scales of Planck length [3. 
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scalar, which are relevant for our work. (We refer the reader to Ref. 3l[ for details). We first note 
that the directional Hubble rates have a universal maxima, given by 

flW = JL . (2.22) 
Using eq. (|2.19p in the definition of the expansion scalar, we obtain 

(sin(/ZiCi +/12C2) +sin(/2 2 C2 +/2 3 c 3 ) +sin(/iici +/2 3 c 3 )) . (2.23) 



2 7 A 



It has a universal maxima given by # ma x = 3/ (27 A). Thus, unlike in the classical GR, where the 
directional Hubble rates and the expansion scalar diverges for arbitrary matter in the Bianchi-I 
spacetime, they are generically bounded in LQC. The boundedness of the directional Hubble rates 
immediately implies the same for the shear scalar a 2 given by eq. (|2.12p . which turns out to be, 



1 



(cos (/x 2 c 2 )(sin (/xici) + sin (/u 3 c 3 )) - cos (/iici)(sin 0« 2 c 2 ) + sin (p 3 c 3 ))) 2 

(2.24) 



3 7 2 A 2 

+ cyclic terms 



with a maxima given by <7 2 max = 10.125/(37 2 A 2 ). Recall that the shear scalar in classical GR 
diverges as 1/V 2 as a singularity is approached. However, in LQC, its behavior is very different in 
the Planck regime where it departs from the classical behavior and is generically bounded without 
any assumption on the matter content. Finally, the expression for energy density, defined as 
p = Hmatt/V, can be obtained from the vanishing of the Hamiltonian constraint, T-L c s ~ 0: 

P = 8n Q 2\2 ( sin (z"i c i) sin (^2C 2 ) + cyclic terms) , (2.25) 
which is also a universally bounded function with a maxima given by 

Thus, in the loop quantization of Bianchi-I spacetime, the expansion and the shear scalar of the 
geodesic congruences and the energy density of the matter content turn out to be generically 
bounded. This is direct consequence of the underlying discreteness of the quantum geometry in 
LQC, captured by A. In the limit, A — > 0, the quantum discreteness disappears and the behavior 
of the above physical quantities turns out to be in agreement with GR, diverging as the volume 
goes to zero. 



III. CURVATURE INVARIANTS 

In GR, one of the characteristic feature of singularities is the associated divergence in the 
curvature invariants. The key question is whether the quantum geometric effects encoded in the 
effective dynamics of LQC lead to a bound on the curvature invariants. The answer to this question 
is not obvious from the properties of energy density, and the expansion and shear scalars, without 
further assumptions on the nature of the considered matter or its equation of state. Since the 
curvature invariants involve the second order time derivatives of the metric, it is possible that even 
when p, 9 and a 2 are bounded, for a certain choice of the equation of state w, curvature invariants 
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may blow up in the evolution. Examples of such events, exist in isotropic LQC, where they originate 
from the divergences in the pressure P of the matter content 33j, |42( . In the following, we study 
the properties of the Ricci scalar R, Kretschmann scalar K where K = R a 0uvR u an d the square 
of the Weyl curvature C^C^^ . n 

In terms of the directional Hubble rates Hi and the second order time derivatives of the direc- 
tional scale factors, the expressions for the Ricci scalar, Kretschmann scalar and the square of the 
Weyl curvature can be written as: 



R = 2\ HtH 2 + H 2 H 3 + H 3 H X + ^ 



i=l 



3 ..o s 

at 



K = A[ HlHl + HfHl + H 2 2 Hl + Y,~ 



(3.1) 



(3.2) 



i=l 



and 



H\H 2 + H 2 H 3 + H 3 H( — H\H 2 H 3 {Hi + H 2 + #3) 



r2 Tj2 



+ — ( H\H 2 — H3H1 + 2H 2 H 3 ) + — ( H\H 2 — H 2 H 3 + 2H 3 H\ 

ai V«i J a 2 \ a 2 



a3 ( a 3 



03 \03 



H 2 H 3 — H 3 H\ + 2H\H 2 



ai _ 02 _ 03 
ai a 2 a 3 



From above expressions, it is clear that since the directional Hubble rates in the effective spacetime 
of Bianchi-I model in LQC are bounded, as evident from eq. (|2.17p . the boundedness of curvature 
invariants is determined by the dynamical equations cij/aj. These can be obtained by using the 
Hamilton's equations for pi (eq. fl2.18j) ) and Cj (eq. (|2.20p ). As an example, the equation for a\/a\ 
can be written as, 



«i 



1 



a\ 27A 



cos(/iiCi - /x 2 c 2 ) ^-(/iici - fi 2 c 2 ) + cosOuxci - fi 3 c 3 ) ^-{mci - /i 3 C 3 ) 
at at 



d 1 _ 5 

+ cos(/i 2 C2 + /U3C3) — (/u 2 c 2 + /U 3 c 3 ) + — - (sin(/xici - p 2 c 2 ) + sin(/xici - ,u 3 c 3 ) + sin(^t 2 C2 + ^303)} 
at 27A 



which for matter with a vanishing anisotropic stress, can be expressed in the following form: 



(3.3) 



°1 = -47rG(p + P) + — ( 6(cos(/xiCi - /i 2 c 2 )K2i + cos(/iiCi - /j, 3 c 3 )k 3 i)H + (k 12 + K13) — 
ai 41/ V pi 



I P2 P3\ \ . 1 



X , 



(3.4) 



11 Other curvature invariants can be studied similarly, or there expressions can be obtained from these invariants. 
As an example, properties of another useful invariant R a pR a ^ can be obtained using the identity: R a pR a ^ = 



9 



where pressure P = 



—dl-l mai tt/dV, and % 2 is a bounded function denned as 



X 



(single: - /x 2 c 2 ) + sin(/iici - /i 3 c 3 ) + sin(/x 2 c 2 + M3C3)) • 



(3.5) 



In order to obtain (|3.4|) . we have used the relation between the connection and triad components 
and the constants Kij (given by eq. (|2.8|) ). Similar calculations, result in the equations for a 2 /a 2 
and 03/03. Using these dynamical equations, along with those for the directional Hubble rates 
(|2.19p . expressions for curvature scalars can be obtained by straightforward calculations. 

Substituting eqs. (|2.17p and (|3.4p . and similar equations for the time derivatives of a 2 and 03 in 
the expression for the Ricci scalar (|3.ip . for matter with a vanishing anisotropic stress, we get, 



R 



-24nG(p + P) + 1 ( ^(«12 + K13) + — («23 + «2l) + — («31 + **2. 
V \Pl P2 P3 



+ 



1 



2 7 2 A 2 



3 + cos 2 (/iiCi) (sin 2 (/i 3 c 3 ) + 4sin(^ 2 c 2 ) sin(^ 3 c 3 ) - cos(2 / u 2 c 2 )) 



+ cos 2 (/i 2 c 2 ) (sin 2 (^ici) + 4sin(/iici)sin(^ 3 c 3 ) - cos(2/x 3 c 3 )) 
+ cos 2 (/i 3 c 3 ) (sin 2 (/i 2 c 2 ) + 4sin(/2ici) sin(/2 2 c 2 ) - cos(2/xici)) 

- (sin 2 (/ZiCi) sin 2 (/i 2 c 2 ) + sin 2 (/ZiCi) sin 2 (^ 3 c 3 ) + sin 2 (^ 2 c 2 ) sin 2 (^ 3 c 3 )) 



(3.6) 



Let us analyze the behavior of Ricci scalar in the effective spacetime. From the dynamical 
equations in Sec. II (eqs. (|2.18p and (|2.25fl ) . p and Pi/pi are universally bounded in LQC. Further, 
are the constants equal to ((HPi — CjPj) in LQC (eq. (12.21j) ). These imply, that the Ricci scalar 
can diverge only if volume vanishes and/or the pressure diverges at a finite energy density. Note 
that if we consider the case of a fluid with a finite equation of state and the dynamical evolution 
is such that the volume does not vanish, then Ricci scalar does not diverge. In such a case, Ricci 
scalar is always bounded. An example of such a matter content is the case of perfect fluids which 
will be discussed in Sec. IVC. 



be 



The computation of the Kretschmann scalar for matter with a vanishing anisotropic stress can 
performed in an analogous way. Using the equations for Hi and 'di/ai, in eq. (|3.2p . we obtain 

x 2 

{(2 + cos(/2ici - P2C2 ~ 2// 3 c 3 ) + cos(/2ici - p 2 c 2 + 2/2 3 c 3 )) sin 2 ( AlC1 ^ 2C2 ) 



K = 

2 7 4 A 4 

+ sin 2 (//ici - P2C2) + 2(cos(/i 2 c 2 ) - cosQiiCi)) sin(/2ici - /2 3 c 3 ) sin(/2 3 c 3 )} 



+ 



1 



4 7 2 A 2 U 2 



— («12 + Kl3j + «23 , 
Pi VP2 



16vrG(p + P)U 



P2 \ _ _ 

+ GH(k 2 i COs(/ilCl - H2C2) + «31 cos(/iiCi - /X3C3)) 

,P2 P3/ 



7A {6# 7 A(k 2 i cos(/2iCi - /2 2 c 2 ) + K31 cos(/2ici - /x 3 c 3 )) 
+ K 23 (cos(/i 2 c 2 )(sin( / uiCi) + sin(/Z 3 c 3 )) - cos(/i 3 c 3 )(sin(^iCi) + sin(/Z 2 c 2 ))) 
+ (/C12 + K13) cos( / uici)(sin(/i 2 c 2 ) + sin(/Z 3 c 3 ))} + 2U (-87rG 7 2 A 2 (p + P) + x 2 ) 
+ cyclic terms . 



(3.7) 



Analysis of various terms in the above expression shows that as in the case of the Ricci scalar, 
the Kretschmann scalar can also diverge in LQC. Potential divergent terms are proportional to 
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1/V 2 , P/V and terms linear and quadratic in pressure. If the dynamical evolution restricts these 
possibilities, then K is bounded. 

Similarly, we can compute the square of the Weyl curvature. In order to obtain its expression, 
it is useful to note that eq. (|3.3p can be written in the following symmetric form 



HtiHx - H 2 ) - H X {{H 2 - H 3 ) 2 + H X H 2 + H 2 H 3 - 2H 2 ] 



+ H 2 (Hf + H 2 H 3 -H 1 H 2 -H 3 H 1 ) + cyclic terms 



(3.8) 



Using the equations for the directional Hubble rates and its time derivative (obtained using (|3.4p ). 
in the above equation, we obtain the following expression for a matter with a vanishing anisotropic 
stress, 



1 



2cos(2(^iCi - /x 2 c 2 ) + cosQuici + fx 2 c 2 ) - cos(2(/xiCi + /i 2 c 2 ))) 



32 7 2 A 2 y 

+ 5 cos(/2iCi - /2 2 c 2 ) + cos(2/2 2 c 2 )(3 - 2 cos(/2 3 c 3 )) + cos(2/i 1 c 1 )(-6 + cos(2/2 3 c 3 )) 
+ 3cos(2/i 3 c 3 ) + 3(sin(2/2ici) — 4cos(/2 2 c 2 )(sin(/2ici) + sin(/2 2 c 2 ))) sin(2/2 3 c 3 ) 

- 2 {8cos 2 (^iCi) sin(/i 2 c 2 ) sin(/x 3 c 3 ) + 2shx{fl\c{) [- cos(2/x 3 c 3 ) sin(/2 2 c 2 ) 
+ 6cos(/z 2 c 2 ) cos(/2 3 C3)(sin(/2ici) + sin(/2 2 c 2 )) - 2 cos 2 (/2 2 c 2 ) sin(/x 3 c 3 )] 

+ 3 cosQuici) (- cos(/2 3 c 3 ) + cos(2/i 2 c 2 + /2 3 c 3 ) + cos(/2 2 c 2 + 2/x 3 c 3 ) 

- 2sin(/2ici)sin( / u 2 c 2 + /x 3 c 3 )) } x — («i 2 + k V3 ) + k 23 f — - — 

J IPl \P2 P3 

+ 6H(k 2 i cos(/iiCi - /x 2 c 2 ) + k 3 i cos(/2iCi - /2 3 c 3 )) - 167rG(p + P)V 



+ 



X 



4 7 4 A 4 



(cos(// 3 c 3 )(sin(//ici) + sin(/i 2 c 2 )) - cos(/uici)(sin(yu 2 c 2 ) + sin(^ 3 c 3 ))) 



(cos(/x 2 c 2 )(sin(^iCi) + sm(// 3 c 3 )) - cos(^iCi)(sin(/z 2 c 2 ) + sin(/i 3 c 3 ))) 



1 



16V < 



-(«3i - «ia) + «23 ( - - - ] + 16vrG(p + P)F 



+ 6H(k 12 COS(^lCl - // 2 C 2 ) - K 3 i COS(/X!Cl - /X3C3)) 



. Pi . P2 \ P3 f 

«12 1 («31 + 

Pi P2 J P3 



- 6H (2«i2 cosmic! - /x 2 c 2 ) - K31 cos(/iiCi - /i 3 c 3 ) - k 23 eos(/i 2 c 2 - ^3)) 
+ cyclic terms . 



(3.9) 



Though more complicated than the expressions for R and K, above expression reveals similar 
properties regarding the boundedness and possible divergences. As for Ricci and Kretschmann 
scalars, the square of the Weyl curvature can also diverge in LQC if the pressure becomes infinite 
and/or volume vanishes at a finite value of the energy density. 



Let us now summarize the main result of this section. Restricting to the case of matter with a 
vanishing anisotropic stress, analysis of the behavior of various curvature invariants reveals some 
important features of loop quantum dynamics in Bianchi-I spacetimes. We have found that in LQC, 
even though energy density and the expansion and shear scalars are universally bounded by values 
determined by eqs, (|2.23)2.24|) and (I2.25p . curvature invariants can potentially diverge. It is to be 
noted that for all matter models, which lead to dynamics such that the pressure is bounded and 
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physical volume never vanishes during the evolution, curvature invariants are bounded. However, 
if there exists a physical solution such that at a finite value of p, 9 and a 2 , pressure diverges and /or 
volume vanishes then curvature scalars can diverge. For conventional matter models, such as dust, 
radiation or stiff matter, these conditions are not satisfied. Conventional singularities for Bianchi- 
I spacetimes for fluids with a finite equation of state, occur with an associated divergence of p, 9 
and a 2 . Since these quantities are generically bounded in LQC, all such singularities are avoided 
in effective spacetime. In all these cases, curvature invariants turn out to be bounded. Thus, in 
contrast to GR, conditions for which curvature invariants can diverge in the effective spacetime 
description of LQC are highly restrictive. This is illustrated in Sec. IVC, where we show that for 
perfect fluids with w > —1, there are no divergences in curvature invariants. On the other hand, 
in GR, such fluids generically lead to divergences in curvature invariant. 

Our result on the curvature invariants generalize the one for the isotropic models in LQC where 
it was demonstrated that curvature scalars can diverge if the pressure becomes infinite at a finite 
value of energy density [jjjj. We find that such singularities can also arise in effective spacetime 
description of Bianchi-I models. Possible divergence in curvature invariants leads to some important 
questions: What is the nature of the singularities associated with events where curvature invariants 
can potentially diverge in Bianchi-I spacetime in LQC? Certainly, these singularities are rather 
special as they occur when p, 9 and a 2 are all finite. However, are these singularities strong or 
weak? And do they imply, breakdown of geodesies? We answer these questions in the next section. 



IV. GEODESICS AND STRENGTH OF SINGULARITIES: GENERAL 
CHARACTERISTICS AND EXAMPLES 

In order to understand the nature of events where curvature invariants diverge, one must analyze 
the behavior of geodesies and the strength of such singular events. It is important to note that 
divergence of curvature invariants though an indication, is not a sufficient condition for a physical 
singularity to exist. As discussed earlier, a singular event in terms of curvature invariants may allow 
a safe passage of detectors, in which case the singularity turns out to be weak. In order to determine 
whether a singularity is strong or weak, we consider the criterion developed by Tipler and Krolak, 
which involves integrals of Ricci and Weyl curvature components over null or time-like geodesies 
[4]. Here, without any loss of generality, we will analyze null geodesies. 12 In the following we first 
analyze the geodesic equations and contrast their behavior in GR and in the effective spacetime 
description of LQC. We then analyze the conditions for the strength of the singularities. Together 
with geodesic equations, these conditions enable us to understand general features of singularity 
resolution in the Bianchi-I spacetime in LQC which point towards the lack of strong singularities. 
We then consider the case of perfect fluids with a vanishing anisotropic stress and a finite equation 
of state w = P/p > —1. Our results show that unlike in GR where such a fluid leads to strong 
singularities of barrel, cigar, pancake and point types, all strong singularities are absent in LQC 
and the effective spacetime is geodesically complete. We also discuss the potential possibilities of 
finite volume singularities which can lead to divergences in curvature invariants. These are possible 
generalizations of sudden singularities in Bianchi-I models 43|, 44 ] . We find that such singularities 



are weak in nature, and geodesies can be extended beyond them. 



Similar analysis can be carried out for time-like geodesies, and similar conclusions as reached in our analysis can 
be obtained. 
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A. Geodesic evolution in Bianchi-I spacetime 

For the metric f)2. 1 1) . the null geodesies obey, 

(« ft )' + r^V = o, (4.i) 

where the 4- velocities u a : u a = dx a /dr, satisfy u a u a = 0, and a 'prime' denotes derivative with 
respect to the affine parameter r. Computing the Christoffel symbols for the metric (|2.1|) . a 
straightforward calculation leads to the following equations: 



x" = -2x't'Hi, y" = -2y't'H 2 , z" = -2z't'H 3 , t" = -a\Hi x' 2 - a 2 2 H 2 y' 2 - a 2 3 H 3 z' 2 (4.2) 
and 

j — ht ,/ - hi ?' — f'-iMx^iM 

of a A 2 a\ \a{ a z 2 a\ 

where ki are constants. From these equations, one can analyze the geodesic extendibility. In 
particular, one finds that the geodesic evolution breaks down at a finite value of the affine 
parameter when any of the directional scale factors Oj — > and/or the directional Hubble rates Hi 
become infinite. 




To analyze the fate of geodesies, we first consider the cases which are expected to arise in 
most physical situations. These correspond to the scenarios where in classical evolution, energy 
density, expansion and shear scalars diverge as the singularity is approached. From eqs. (|3.H - 
3.3j) . we see that such divergences cause curvature invariants to blow up. These cases have been 
extensively studied for the perfect fluids [49-54]. In terms of the directional scale factors (and their 
permutations), the resulting singularities are of the following types: (i) barrel, characterized by 
a\ — > finite value, 02,03 — > 0, (ii) cigar, where a\ — > 00, 02,03 — > 0, (iii) pancake, where ai — > 0, 
a 2j03 — > finite value, and (iv) point or isotropic, characterized by all scale factors vanishing. For 
all these singularities, at least one of the scale factors and the physical volume vanishes at the 
singularity. Since the directional Hubble rates diverge at these singularities in the classical theory, 
eqs. (|2.11U2.12l[2TT4"|) and eqs. (|3.HI3.3p . imply the divergence of p, 9, a 2 and the curvature invariants 
in the classical theory. From eqs. (|4.2p and (|4.3p . we find that the geodesic equations break down 
at above singularities in the classical theory. 

On the other hand, in LQC, for all such cases geodesic equations do not break down. Classical 
singularities which occur such that Oj — > and \Hi\ — > 00 simultaneously, are forbidden in LQC. 
Since Hi is bounded (eq. (|2.22p ). the scale factors at which above singularities occur are excluded 
from the effective spacetime. To illustrate, in a typical evolution where a resolution of above 
singularities occurs, one starts with an initial data which for the classical evolution leads to a 
vanishing of a scale factor and a divergence of the corresponding directional Hubble rate in a 
finite value of affine parameter. For concreteness, let these be a\ and Hi respectively, where 
we assume Hi to be positive and increasing in the backward evolution towards the classical 
singularity. Let us first note, that when the spacetime curvature is small (fliCi ^ 1), effective 
dynamical equations approximate the classical dynamical equations. This can be explicitly seen 
by comparing the effective Hamiltonian constraint (12.150 with the classical Hamiltonian constraint 
(|2.4p . In this regime, Hi in LQC approximates its counterpart in GR. However, as Hi becomes 
large, spacetime curvature increases, and deviations between LQC and GR become significant. 
Analysis of eq. (|2.17p and (|3.4p shows that depending on the values of fiiCi in the evolution, Hi 
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attains a maximum value in the Planck regime, and starts decreasing, subsequently vanishing 
before becoming negative. 13 This causes a turn around of the scale factor a\ in effective dynamics 
and the singularity is avoided. 14 Unlike GR, where H\ diverges and a\ approaches zero, effective 
dynamical equations cause a turn around of a\ and the vanishing of a\ does not occur. Similar 
conclusion holds, if more than one scale factors vanish at a classical singularity with associated 
divergences in the directional Hubble rates. Thus, due to quantum gravitational modifications 
encoded in effective dynamics, classical singularities accompanied by a% — > and a divergence in 
Hi are avoided. Eqs. (|4.2p and (|4.3|) . then immediately imply that geodesic evolution in these cases 
does not break down in the effective spacetime description of Bianchi-I spacetime in LQC. 



Remark: In GR, divergence of p, Hi, 9 and a 2 may not always occur at a vanishing physical vol- 
ume. In isotropic models, singularities where p and isotropic Hubble rate diverge at non-vanishing 



volume have been investigated recently (see for eg. 35|, [4l|. These are the big rip (occurring at 
infinite volume) and big freeze (occurring at a finite non-zero volume) singularities, which are 
geodesically inextendible events in GR. In isotropic spacetimes in LQC, such singularities have 
been shown to be avoided due to the boundedness of Hubble rate in LQC 0,0, 0,0- In the 
effective dynamics, Hubble rate vanishes in the Planck regime and causes a turn-around of the 
scale factor. This leads to a recollapse of the universe before scale factors at which above classical 
singularities occur. Assuming that such singularities exist in Bianchi-I models in GR, these will 
result in break down of geodesies in the classical theory due to divergences in the directional Hubble 
rates. Since Hi are bounded in LQC, above singularities will be avoided in LQC in the same way 
as singularities at vanishing scale factors (with associated divergence in Hi) are avoided. From 
eqs. (|4.2p and (|4.3p . we find that the geodesic equations will be well behaved for these cases in LQC. 



So far we have discussed the cases of singularities where divergence in curvature scalars is 
accompanied by a divergence in p, 9 and a 2 . For all such cases, we have found that geodesic 
evolution does not break down in effective spacetime description of LQC. We now turn to the 
cases where curvature invariants can diverge at a finite value of p, 9 and a 1 . Examples of such 
singularities are not known in Bianchi-I models. 15 Here we would assume that such possibilities 
can arise in the physical evolution in LQC. From the previous section, we recall that curvature 
invariants in LQC can diverge under following conditions: if the physical volume of the spacetime 
approaches zero, and/or the pressure becomes infinite in magnitude. If these conditions are 
satisfied for physical solutions of effective dynamics in LQC then following two contrasting 
scenarios can arise: 

(i) If in the physical evolution, curvature invariants diverge at a finite value of p, 9 and a 2 , 
and a vanishing volume, (which implies that at least one of the scale factors vanishes) at a 
finite value of the affine parameter, then geodesic equations (|4.2j) and (|4.3p break down. If such 
curvature invariant diverging events exist, then geodesic evolution is not complete. 16 

(ii) If in the physical evolution, pressure becomes infinite at a finite value of p, 9 and a 2 , and 



13 For a pictorial illustration of this behavior in a model in effective dynamics in Bianchi-1 spacetime, see Ref. [27( 
(Fig. 4). 

14 The situation is analogous to isotropic models in LQC, where in the backward evolution towards the big bang, a 
bounce of the scale factor occurs after the Hubble rate attains its maximum value. For examples, see Ref. [34l. log]. 

15 In the anisotropic models in classical theory, such singularities occurring at finite scale factors have been investi- 
gated in Bianchi-VIIo spacetime [44| . 

16 However, there may still be no physical singularity. An example being if the spacetime is maximally extendible. 
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a non-vanishing values of the scale factors, and hence non-vanishing volume, then the geodesic 
equations are well behaved. For these curvature invariant diverging events, geodesic evolution is 
complete in the effective spacetime description of LQC. The behavior in this case is similar to the 



case of pressure singularities occurring in isotropic models in LQC [33|, |42|, |48|], where though for 

an exotic choice of matter, curvature invariants diverge in the evolution, but the spacetime turns 
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out to be geodesically complete [33[] 



To summarize this part, we find that irrespective of the choice of matter, geodesic evolution in 
the effective spacetime description of Bianchi-I LQC is complete for all the cases except the isolated 
case when physical evolution may allow divergence of curvature invariants at finite values of p, 9 
and a 2 with at least one of the scale factors vanishing, at a finite value of the affine parameter. No 
known classical singularities of Bianchi-I model satisfy the conditions for this isolated case, and at 
this stage, its existence is only a potential possibility. 



o 



B. Strength of singularities 

The strength of singularities can be determined using the necessary and sufficient conditions 
obtained by Tipler [1] and Krolak These conditions are used to classify the singularities as 
strong and weak types, which provide insights on the magnitude of tidal forces experienced by an 
in-falling detector (or an observer) towards the singularity 0]. According to the Krolak's criteria, 
the conditions to determine the strength of the singularities are the following. If for a null geodesic, 
the integral over Ricci curvature components, 

dr'R^ u^u u (4.4) 

or the integral over Weyl curvature components 

dT"\C a0flv u^u v \\ (4.5) 

diverges, then the singularity is considered to be strong. Else the singularity is weak. Tipler's 
conditions are similar to those of Krolak, but involve an additional integral over the affine 
parameter Q]. As an illustration of when a singularity is strong or weak, if in a physical evolution 
R a pu a u^ is proportional to l/r m and C^^u^u 11 is proportional to l/r n , then a singularity is 
strong curvature type by Krolak's criteria if m > 1 or n > 3/2. Similarly, a singularity is strong 
following Tipler's criteria if m or n > 2. Thus, a singularity can be strong by Krolak's criteria, 
yet it can be weak according to Tipler's criteria. However, all singularities which are strong by 
Tipler's criteria are also strong by Krolak's criteria. Due to this reason, we will consider Krolak's 
criteria to address the resolution of strong curvature singularities in our analysis. 

To analyze the strength of potentially singular events in effective spacetime description of 
LQC, we find the expressions for the integrands in the integrals (|4.4|) and (I4.5p . For the met- 
ric (|2.ip . the only non- vanishing components of Ricci and Weyl tensors are Ra (i = 1..4) and 
Ci2i2> C1313, Ci4i4, C2323, C2424, C3434 (and those obtained by symmetric and anti-symmetric trans- 
formations on the indices, such as (^2121^2112 etc.). Using eq. (|4.3p with the expression for the 




17 Geodesies can be extended beyond such events even in GR [45 
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Ricci tensor components in terms of first and second order time derivatives of the directional scale 
factors di, one can obtain the integrand for (|4.4p as 



u u 



-f ( H1H2 + H1H3 ) + -| ( H1H2 + H2H3 



A- 



m a 3 



a.2 «3 

H — I ( H 3 H 1 + H 2 H 3 — -i — "— 
a| V ai a 2 

where fcj are constants obtained from the geodesic equation f|4.3[) . 
A similar calculation for the integrand of (|4.5p yields, 



ai a 3 



Caftcd U a U b 



6a? 



(1 + a?) ( #i# 2 + #i# 3 - 2H 2 H 3 + — + — - 2— ) + 3(a| - a| 

a 2 a 3 oW " ' V a 3 



«3 



h H1H2 — H\H 3 

a-2 



a.3 „d2 



+ 



k 2 
K y 

6a| 



+ — - 2— ] + 3(a^ - af) tf 2 F 3 - #i# 2 + 

a 3 a 2 J \ ai a 3 



(1 + ( H 2 H 3 + H X H 2 - 2H X H 3 + ^ 



ai a 3 



+ — 

«3 



(1 + a z 3 ) ( # 2 tf 3 + #i#3 - 2H 1 H 2 + — + — - 2— 

ai a 2 03 



2 2 \( fl l a 2 
+ 3(a 2 - ai) 



+ H2H3 — H\H 3 



ai a 2 

11 1 ( ky k\ 
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x I 2 2 "t" 2 

«I «2 «3 

,2 1.2 
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(4.6) 
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(4.7) 



In order to determine whether a singularity is strong in effective spacetime description of LQC, 
we substitute (|4,6p and (|4.7|) in (|4.4p and (|4.5p respectively. Both of the above integrands, involve 
quadratic terms in the directional Hubble rates, linear terms in di/ai and terms with inverse scale 
factors. Expressing Oj/a, in terms of Hi and Hi, integration by parts of (14. 4|) and (|4.5p yields 
integrals with terms with directional Hubble rates and 1/a?. Since, directional Hubble rates are 
bounded functions (I2.22p . whether or not the integrals in Krolak's conditions diverge depend on 
the behavior of 1/a? terms. In particular, if in an evolution a, — > as the singularity is approached 
in a finite value of affine parameter, then integrals (|4.4h and (14. 5p can diverge, leading to a strong 
singularity a la Krolak. 



Let us determine the strength of those curvature invariant diverging events in LQC which occur 
due to divergence in pressure at finite non-vanishing value of directional scale factors (or finite 
volume). For such 'pressure' singularities, divergence of curvature invariants appears only due to 
divergence in di/ai terms. Since these terms are integrated over affine parameter at least once in the 



16 



Krolak's conditions (and at least twice in the Tipler's) conditions, they do not yield divergences of 
the integrals (|4.4p and (|4.5p . Pressure singularities occurring at a finite value of scale factors, thus, 
turn out to be weak singularities in the effective spacetime description of Bianchi-I model in LQC. 
This conclusion supplements the result obtained in the previous sub-section where we established 
that for such singularities, geodesic evolution does not break down. Thus, these singularities are 
harmless - they are weak in strength and do not forbid extension of geodesies. We recall that 
the properties of pressure singularities occurring at finite scale factors are similar to the isotropic 
models in LQC, where as for Bianchi-I model, they turn out to be weak and geodesically extendible 




The only case where strength of the singular events in LQC can be strong is when at least one of 
the scale factors vanishes at a finite value of p, 9 and a 2 in the physical evolution at a finite value of 
affine parameter. This is the identical condition for the isolated case of geodesic inextendibility for a 
mathematically allowed possibility in effective dynamics of LQC found in Sec. IVA. If such an event 
turns out to be physically realizable in effective dynamics of LQC, then one can explicitly compute 
the integrals (|4.4j) and (14. 5h and determine whether it corresponds to a strong or a weak singularity. 

In conclusion, we find that those events which lead to a divergence in curvature invariants in 
LQC are weak singularities if they occur due to divergence in pressure at non-vanishing scale factors. 
If the divergence in curvature invariants is associated with vanishing of one more scale factors at 
a finite energy density and expansion and shear scalars at a finite value of affine parameter, then 
the singularity can be strong or weak. We emphasize that examples of events of latter type in GR 
or LQC are not known, and this case may turn out to be physically not realizable. 



C. Fate of geodesies and strong singularities: physical examples 

We now discuss specific examples where the implications of modified loop quantum dynamics 
on singularity resolution can be manifestly seen. We start with the case of a perfect fluid with an 
the equation of state w > — 1 and a vanishing anisotropic stress. Perfect fluids with different values 
of w have been extensively studied in the Bianchi-I models in GR and in many situations analytical 
solutions are also known 4^-54]. Singularities in the classical theory, which are of barrel, cigar, 



pancake and point like, are strong and geodesically inextendible events occurring at a vanishing 
volume. We will show that in LQC, all these singularities are resolved. We then discuss aspects of 
singularity resolution for potential curvature invariant divergent events in Bianchi-I model in LQC 
which can arise for more general fluids satisfying a generalized equation of state P = P(p), which 
in principle allow singularities other than those found for perfect fluids in GR. These would be 
anisotropic generalization of big rip, big freeze and sudden singularities found in isotropic models 
*35l . 41, 43]. Since explicit examples of equations of state which lead to these singularities in Bianchi- 



I model are not yet known, our discussion will only point out the general features of resolution of 
such potential singularities. 



1. Perfect fluid with w > — 1 

A perfect fluid with a vanishing anisotropic stress and a constant equation of state w = P/p 
which is greater than -1 constitutes a large type of matter models which have been studied in 
Bianchi-I spacetimes in GR, which include dust (w = 0), radiation (w = 1/3) and stiff matter (w = 
1) 0-[54]. In the effective dynamics of Bianchi-I model in LQC, quantum geometric modifications 
do not influence the matter part of the Hamiltonian constraint, and the stress-energy tensor of the 
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matter content satisfies the conservation law T%;^ = 0, where 

T„ v = (P + P)v„v u + Pg^ . (4.8) 

The conservation law leads to 

p+{H 1 + H 2 + H 3 )(l + w)p = , (4.9) 



which on integration gives, 



p = C(aia 2 a 3 r (1+w) (4.10) 



where C is a constant determined by the initial conditions. 

Let us recall some features of the classical evolution for the perfect fluids. For w > — 1, energy 
density diverges as the mean scale factor a = (oj^as) 1 / 3 approaches zero in a finite proper time. 
At these events directional Hubble rates diverge, causing expansion and the shear scalar to become 
infinite in GR. Due to these divergences, expressions for Ricci (|3.1|) . Kretschmann ([3.2p and the 
square of the Weyl curvature (|3.3p reveal that these curvature invariants grow unbounded as the 
singularity in classical theory is reached. The singularity, depending on w and initial anisotropics 
can be of the form of a barrel, cigar, pancake or a point [49h51|]. Also, integrals (|4.4l I4.5|) to 



determine the strength of singularities diverge and the singularities turn out to be strong. From 
the geodesic equations (|4.2p and (|4.3p . we find that at these events, geodesic evolution breaks 
down. Thus, in classical GR, perfect fluids with w > — 1 generically lead to physical singularities 
in Bianchi-I model. 

We now analyze the existence of above classical singularities in LQC. Since the conservation law 
(|4,9p is unmodified in LQC, the resulting proportionality of energy density with the scale factors, 
turns out to be identical as in GR. However, unlike GR, where p can grow unboundedly as V — > 
for w > —1, it is bounded by a universal value in LQC. Using the upper bound on energy density 
fl£Zg) , eq.flHU]) yields, 

^«>«3)-< 1+ «> < — ^ ■ (4.11) 

This inequality implies that in a physical evolution, volume V = 010203 never becomes smaller 
than a minimum value given by, 

Fmi „ = p^)*. (4 , 2) 

Given that P = wp, the bound on energy density also leads to a bound on the pressure of the 
perfect fluid: P max = wp m£LX 0Alwpp\. Recall that in Sec. Ill we showed that the curvature 
invariants in LQC can diverge only if the physical volume vanishes or the pressure of the matter 
content diverges. Using above results in eqs, (|3.6l I3.7P and (|3.9p . we find that for a perfect fluid 
with a finite equation of state w > — 1, Ricci and Kretschmann scalars and the square of the Weyl 
curvature are always bounded in LQC. 

The existence of a minimum volume (|4.12p in LQC, leads to another important implication. 
Note that all the classical singularities for a perfect fluid with w > — 1, whether they are barrel, 
cigar, pancake or point like, occur at a vanishing physical volume 43-54]. Irrespective of the 



geometry of the singularity, a non- vanishing physical volume implies that the values of the scale 
factors where these singularities occur in GR are excluded by the effective dynamics in LQC. 
Along with the bounds on the directional Hubble rates, Hi max = 3/(2 7 A) obtained from ([2T7j) . 
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this ensures that the geodesic equations f|4.2j) Sz (|4.3p remain well defined in LQC. Thus, for a 
perfect fluid with w > — 1 geodesic evolution never breaks down in LQC. This result is in sharp 
contrast to the behavior of geodesies for perfect fluids with w > — 1 in GR, where irrespective of 
the choice of w or initial conditions, classical spacetime is geodesically incomplete. 

Finally, let us consider the integrands (|4.6p and (|4.7p . Since the vanishing of the scale factors 
a,i in the effective dynamics of LQC is ruled out for perfect fluid with w > — 1 and the directional 
Hubble rates are bounded, these integrands can only diverge if di/ai become infinite. However, 
eqs. (|3.4p (and similar equations for 02/02 and 03/03) imply that for a non- vanishing volume and 
finite pressure, di/ai are always bounded in LQC. Therefore, the integrands (14. 6p and (14. 7p are 
always bounded. The integrals f|4.4j) and f|4.5|) over these integrands are thus finite for any finite 
range of integration. Hence Krolak and Tipler's conditions for the existence of strong singularities 
fail to be satisfied. We thus conclude that for a perfect fluid with w > —1, geodesic evolution is 
well defined for all times in LQC and there is neither a divergence of curvature invariants nor any 
strong singularities in the effective spacetime. 

2. Potential exotic singularities 

For a perfect fluid with w > — 1, singularities in classical theory occur when physical volume 
vanish. However, if we consider a fluid with a generalized equation of state, singularities at 
finite volume are potentially possible. 18 These singularities come in variety of forms and have 
been classified depending on the properties of Hubble rate, d/a and higher order time derivatives 



67h69||. Since the Krolak's and Tipler's conditions require at least the second time derivative 
of scale factor to diverge in order for a singularity to be potentially strong, we will restrict our 
discussion to only such events, i.e. to big rip, big freeze and sudden singularities. In the following, 
we assume that generalization of these singularities exist in the Bianchi-I spacetime in GR and 
discuss the effects of loop quantization on them. As in the previous section, we assume a fluid 
with a vanishing anisotropic stress. 

(i) Big rip singularities: In isotropic models in GR, this singularity occurs when the scale 
factor diverges with a divergence in energy density and Hubble rate. In LQC, for isotropic 



models, these singularities have been shown to be absent [33l . 1361 . 1371 . I42l |. Its generalization in 
Bianchi-I spacetime would involve a divergence of at least one of the scale factors, directional 
Hubble rate and energy density, causing 9 and a 1 to become infinite. Due to these divergences, 
integrals (|4.4p and (|4,5p would diverge, and these singularities will be strong curvature type 
in GR (a la Krolak's criteria). Analysis of geodesic equations (|4.2|4.3|> shows that they would 
break down in the classical theory, due to divergence in directional Hubble rates. On the 
other hand, in the effective dynamics of LQC, energy density and directional Hubble rates are 
bounded, and hence the scale factors at which these divergences occur will be excluded from 
the effective spacetime. Thus, big rip singularities of Bianchi-I models in GR will be absent in LQC. 

(ii) Sudden singularities: These singularities occur at a finite value of scale factor, energy density 
and Hubble rate but a divergent pressure in isotropic models (43|. In the anisotropic setting, a 
sudden singularity would occur under similar conditions: a divergence in pressure at a finite value 
of the mean scale factor and energy density. From eq. (|3.4p we find that this divergence would cause 
di/ai to blow up, resulting in curvature invariants (j3.6H3.9p to become infinite. Since these singu- 



18 For examples of models, in isotropic spacetimes in GR, where such generalized equations of state have been 
extensively studied, we refer the reader to Ref. [35I. l4lL 143. l67l - [69t . 
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larities would occur at finite cjj > and with a divergence only in the pressure, analysis of Sec. IVB 
shows that these will be weak singularities, both in GR as well as LQC. Analysis of geodesic equa- 
tions (|4.2j) and (|4,3p shows that as in the isotropic case, these events will be geodesically extendible. 

(hi) Big freeze singularities: These singularities share the properties of big rip and sudden singu- 
larities. In isotropic models, these occur at a finite scale factor but with a divergence in energy 
density and Hubble rate. They are strong and geodesically inextendible in GR. In LQC, such 
singularities have been shown to be forbidden in isotropic models 33|, 42]. Their fate in Bianchi-I 



models will be similar. Due to bounds on directional Hubble rates and expansion scalar (|2.17p and 
(|2.23p . effective dynamics would cause a turn-around of scale factors before the classical singularity 
could be reached, and these singularities would be avoided. 

In summary, if we allow generalizations of big rip, sudden and big freeze singularities in Bianchi- 
I spacetime in GR, then of these only sudden singularities are the ones that are not excluded in the 
effective spacetime description of LQC. Though at these events, curvature invariants in effective 
spacetime description of LQC diverge, these turn be weak singularities beyond which geodesies can 
be extended. 



V. SUMMARY 

It is widely believed that a viable theory of quantum gravity must address the way space-like 
singularities of the GR are overcome. An important issue is whether such a theory resolves all 
space-like singularities, in particular those beyond which geodesies in classical theory can not be 
extended, and which are strong by strength. 19 The goal of this work, is to take a first step in 
answering this question for the loop quantization of Bianchi-I model. In recent years, resolution 
of big bang/crunch singularities has been successfully demonstrated in LQC for different models, 
and detailed analytical, phenomenological and numerical aspects of the quantum theory of various 
cosmological spacetimes have been studied. An effective spacetime description of the underl yin g 



quantum theory has also been derived, usin g ge ometrical formulation of quantum mechanics 38] 



It is based on an effective Hamiltonian 6114631]. whose resulting dynamics turns out to have an 
excellent agreement with the quantum evolution in isotropic and anisotropic models [lMIGIGI 
40J]. Effective dynamics has been extensively used to extract various novel physical predictions in 
LQC (see Sec. V of Ref. [8'] for a review). These developments in LQC, set the stage to investigate 
the issue of generic resolution of strong singularities in Bianchi-I models. 

To understand the general nature of singularity resolution, key properties of curvature invari- 
ants, geodesic equations and strength of any singular events need to be analyzed. These were 
investigated in this work, assuming the validity of effective spacetime description, and for matter 
with a vanishing anisotropic stress. Using effective dynamics of Bianchi-I model, it has been earlier 
shown that for generic matter, energy density (p) and directional Hubble rates (Hi) are bounded 
by universal values. The latter lead to bounds on expansion (9) and shear (a 2 ) scalars of geodesic 



congruences in the effective spacetime [23, |31|, |32| • These bounds are direct ramifications of the 
underlying quantum geometry in LQC. Our results on the properties of curvature invariants show 
that though quantum geometric effects bind them in many cases of physical interest, interestingly, 
they can in principle diverge. Thus, the bounds on p, 9 and a 2 do not guarantee the boundedness 
of curvature invariants. We show that curvature invariants are bounded in Bianchi-I model in 



Questions about the resolution of other singularities can also be raised in quantum gravity. As an example, it has 
been argued that a viable theory of quantum gravity should not resolve a certain class of time-like singularities 
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LQC except if in the physical evolution, pressure becomes infinite and/or the volume vanishes. 
Since these conditions must be satisfied at a finite value of p, 6 and a 2 , divergences of curvature 
invariants occur under very special conditions. In isotropic models, for a generalized equation of 
state, divergence in pressure at finite scale factor and energy density results in a sudden singularity 
which is weak and a geodesically extendible event in classical theory as well as LQC [33j]. General- 
ization of these singularities in Bianchi-I models, satisfy the conditions for divergence in curvature 
invariants in LQC. However, unlike pressure singularities, examples of potential singularities which 
may occur at vanishing volume with a finite energy density and expansion and shear scalars are 
not known. 

Existence of events where curvature invariants can potentially diverge, does not necessarily im- 
ply existence of a strong singularity or geodesic inextendibility. We show that this turns out to 
be true in the Bianchi-I model in LQC. Analysis of the (null) geodesic equations in the effective 
spacetime reveals that geodesic evolution remains well defined for events where curvature invari- 
ants diverge due to pressure becoming infinite at finite scale factors. Such singularities also turn 
out to be weak, and thus are harmless. These curvature invariant divergent events, thus do not 
lead to physical singularities in LQC. Only when the divergence in curvature invariants occurs 
at a vanishing volume with a finite value of energy density, and expansion and shear scalars, at 
a finite value of affine parameter, can the geodesic equations in effective spacetime of Bianchi-I 
model in LQC break down. None of the known classical singularities of Bianchi-I models, satisfies 
this condition. At this stage, this isolated case is only a potential possibility which may not be 
physically realized in the effective dynamics. Assuming that such a curvature invariant diverging 
event is allowed by physical solutions, the resulting singularity may be strong or weak, depending 
on the details of the dynamical evolution. To understand the existence of such a potential sin- 
gularity, it will be important to understand the validity of effective spacetime description when 
scale factors approach zero. If effective dynamics gets additional corrections at small scale factors, 
above potentially curvature invariant diverging event may be ruled out. 20 It is important to note 
that such a singularity does not arise for matter which obeys positive energy conditions. Hence, it 
is also possible that such a potential singularity is realized for unphysical conditions for matter in 
the effective dynamics, which may be restricted by the quantum theory. 

These results stand in sharp distinction to the ones in the classical theory where Bianchi-I 
models generically lead to strong singularities and the spacetime is geodesically incomplete. To 
illustrate this, we consider a perfect fluid with equation of state w > — 1. In classical theory, this 
choice leads to barrel, cigar, pancake and point singularities which are strong and geodesically 
inextendible In contrast, we show that all these singularities are resolved in LQC. 



Curvature invariants turn out to be bounded and geodesic evolution never breaks down for perfect 
fluid with w > — 1 in LQC. We also discussed the potential generalizations of big rip, sudden and big 
freeze singularities in Bianchi-I spacetime. Analysis of the general features of these singularities 
using effective dynamics of Bianchi-I model shows that big rip and big freeze singularities are 
resolved by quantum geometric effects, whereas sudden singularities turn out to be of weak strength. 
As discussed above, these will correspond to pressure singularities in Bianchi-I LQC where curvature 
invariants diverge. However, geodesies are extendible beyond them. 

Analysis in this work generalizes the one for the spatially flat isotropic models in LQC performed 
earlier [33|. There we found that all strong singularities were resolved and the effective spacetime 
turned out to be geodesically complete. Analysis of the effective dynamics for k = ±1 model, 



suggests the same result for spatially curved models 42J. We have now shown that similar results 



20 If the topology is considered compact, terms originating from inverse scale factor corrections in loop quantization 
(see Ref. @] for details), would also potentially play an important role in this regime. 
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appear in the presence of non-vanishing Weyl curvature. To understand the physical implications 
of our analysis in more detail, it is important to include the cosmological constant and matter with 
equation of state w < — 1. It is also important to study the implications for the vacuum Bianchi-I 
model where the singularities result only from the divergence in Weyl curvature. Though these 
can be considered as a straightforward generalization of our analysis, given the rich structure of 
singularities in classical Bianchi-I model, important subtleties can not be ruled out. Further, our 
analysis assumed matter with a vanishing anisotropic stress. In future work, this assumption will 
be relaxed. Though we do not expect a qualitative change of results, presence of matter with a 
non-vanishing anisotropic stress is expected to enrich the phenomenological implications. 

The genericity of these results on singularity resolution, achieved earlier for isotropic models 
[33L l42| and here for Bianchi-I model, suggest existence of a non-singularity theorem in quantum 
gravity analogous to the singularity theorems in classical theory. Though these works can regarded 
as a first steps towards realizing such a non-singularity theorem, in future works, several steps need 
to be carefully and systematically taken as we consider more complex models to understand the 
underlying conditions of such a theorem. Some insights on these conditions arise from the recent 



work on detailed contrasts on physics of anisotropic models performed in Ref. 31] . There it was 



found that inclusion of spatial curvature in anisotropic setting, as in the case of quantization of 



Bianchi-II and Bianchi-IX models done in Refs. 23|, [24|, reveals two interesting features. First of 



these is the non-trivial role of energy conditions, and second is the role of inverse triad (or volume) 
corrections (in spatially curved models) on the physics of singularity resolution. Both of these 
features are expected to play an important role in proving such a theorem. These observations also 
make our expectations on genericness of non-singularity results stronger. First, it is well understood 
in LQC (see for eg. Ref. [71[), that inverse volume modifications for model with positive curvature 
help the singularity resolution. 21 Thus, we expect that with modifications coming from both the 
holonomies yielding trignometric functions in the Hamiltonian constraint and the inverse triad 
effects, singularity resolution results would strengthen. Second, even in GR, singularity theorems 
are proved by demanding that matter satisfies certain energy conditions. Energy conditions can 
weed out unphysical solutions, such as those which violate weak energy conditions. Thus, inclusion 
of energy conditions which restrict unphysical matter is expected to bring us closer to prove a 
non-singularity theorem on similar lines as the singularity theorems in the classical theory. 

Next steps in the above direction to prove genericity of singularity resolution would require 
going beyond the effective spacetime description, and the inclusion of inhomogeneities. For the 
first step, one will require a derivation of these results in full loop quantization. Though techniques 
exist to address these issues at the quantum theory level for some matter, these would require 
extensions to arbitrary matter or with certain energy conditions (such as the weak energy condi- 
tion) to prove genericity of singularity resolution. So far little work has been done on the inclusion 
of inhomogeneities in loop quantization of symmetric models, though useful insights have been 



gained by treating inhomogeneities a la Fock quantization methods 64]. Such a hybrid method, 
suggests resolution of singularities in a qualitatively similar fashion as in isotropic and anisotropic 
spacetimes in LQC. In future work, it will be important to understand various aspects of singu- 
larity resolution in such hybrid methods for arbitrary matter, both at the effective spacetime and 
quantum theory levels. Perhaps the most challenging step would then be to understand the full 
loop quantization of inhomogeneous spacetimes. In this direction, important insights are expected 



to arise by carrying forward works attempting to link LQG with LQC 73|, |74|, and also by under- 



standing detailed physics of incorporation of inhomogenities in LQG using spinfoam methods 75] ■ 



21 In fact, it is possible to obtain bounce in As = 1 model with inverse volume modifications, even if holonomy 
modifications are not considered (see for eg. [72 
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Given the nascent stage of results on inclusion of inhomogenities in LQC and relation with LQG, 
at this stage, one can only speculate about the possible roadblocks on the road to prove a non- 
singularity theorem. On general grounds, we expect one of these to be related with the difficulties 
in expressing holonomies as almost periodic functions of connection, which can make analysis of 
the physics of singularity resolution more involved. Fortunately, quantization of Bianchi-II 23] and 



Bianchi-IX spacetimes 24 1, performed without the feature of almost periodicity, and the resulting 

311 ]. provide important hints to proceed in such a situation. It is hoped that 



analysis of new physics 

future work on symmetry reduced models and on inclusion of inhomogeneities, on the lines of the 
analysis performed here, will provide vital clues and a deeper understanding on this fundamental 
issue in quantum gravity. 
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